Astronomy & Astrophysics manuscript no. G3520'astroph 


©ESQ 2013 


March 19,2013 





Letter to the Editor 



m 
o 

(N 



oo 

< 
o 

6 

CO 



> 

(N 

rn 
O 
m 



X 



A candidate circumbinary Keplerian disk in G35.20-0.74N: 

A study with ALMA 

A. Sanchez-Monge 1 , R. Cesaroni 1 , M. T. Beltran 1 , M. S. N. Kumar 2 , T. Stanke 3 , H. Zinnecker 4 , S. Etoka 5 - 6 , D. Galli 1 , 
C. A. Hummel 3 , L. Moscadelli 1 , T. Preibisch 7 , T. Ratzka 7 , F. F. S. van der Tak 8 - 9 , S. Vig 10 , C. M. Walmsley 1 - 11 , and 

K.-S. Wang 12 



INAF, Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 1-50125 Firenze, Italy 
Centro de Astroffsica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 
ESO, Karl-Schwarzschild-Strasse 2, D-85748 Garching bei Miinchen, Germany 

SOFIA Science Center, NASA Ames Research Center, Mailstop 232-12, Moffett Field, CA 94035, USA 

Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK 

Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany 

Universitats-Sternwarte Miinchen, Ludwig-Maximilians-Universitat, Scheinerstrasse 1, 81679 Miinchen, Germany 
SRON Netherlands Institute for Space Research, P.O. Box 800, 9700 AV, Groningen, The Netherlands 
Kapteyn Astronomical Institute, University of Groningen, The Netherlands 

Department of Earth and Space Science, Indian Institute of Space Science and Technology, Thiruvananthapuram, India 
Dublin Institute for Advanced Studies (DIAS), 31 Fitzwilliam Place, Dublin 2, Ireland 
Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 



Received date; accepted date 



ABSTRACT 



We report on ALMA observations of continuum and molecular line emission with 074 resolution towards the high-mass star forming 
region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In 
one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of ~18 M Q . This 
finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3 fim towards G35. 20-0.74 N. 
The disk radius and mass are ^2500 au and ~3 M Q . To reconcile the observed bolometric luminosity (~3 x 10 4 L Q ) with the estimated 
stellar mass of 18 M Q , we propose that the latter is the total mass of a binary system. 

Key words. Stars: formation - ISM: individual objects: G35.20-0.74 N - ISM: jets and outflows 



1. Introduction 

While different scenarios have been proposed to explain the for- 
mation of high-mass (i.e. OB-type) stars (monolithic collapse in 
a turbulence-dominated core - Krumholz et al. 2009] competi- 
tive accretion driven by a stellar cluster - Bonnell & Bate 2006; 
Bondi-Hoyle accretion - Keto I2007I see also the review by 
Zinnecker & Yorke 2007), all of them predict the formation of 
circumstellar disks. It is thus surprising that only a handful of 
disk candidates have been observed in association with massive 
(proto)stars. As a matter of fact, despite many observational ef- 
forts, convincing evidence of disks has been found only around 
early B-type (proto)stars, while circumstellar disks around O- 
type stars remain elusive (Wang et al. I20121 Cesaroni et al. I2007I 
and references therein). Moreover, a detailed investigation of the 
disk properties, comparable to that performed in disks around 
low-mass stars (e.g. Dutrey et al. 2007) is still missing due to 
the large distances of OB-type (proto)stars and the limited an- 
gular resolution at (sub)millimeter wavelengths. With the advent 
of the Atacama Large Millimeter Array (ALMA) the situation is 
bound to improve dramatically, as resolutions <k1" will be easily 
obtained. 



With this in mind, we performed ALMA Cycle observa- 
tions of two IR sources containing B-type (proto)stars. These 
were chosen on the basis of their luminosities (on the order of 
10 4 L Q ), presence of bipolar nebulosities/outflows, detection of 
broad line wings in typical jet/outflow tracers (SiO), and strong 
emission in hot molecular core (HMC) tracers (such as methyl 
cyanide, CH3CN). Here we present the most important results 
obtained for one of the two sources, G35.20-0.74N. 

G35.20-0.74N is a well known star forming region located 
at a distance of 2.19+^1 kpc (Zhang et al. l20Tj9i l. with a lu- 
minosity of ~3 x 10 4 LcB The region is characterized by the 
presence of a butterfly shaped reflection nebula oriented NE- 
SW (see Fig. [T), as well as a bipolar molecular outflow in the 
same direction, observed in 12 CO by Dent et al. ( 1985a), Gibb et 
al. (120031 hereafter GHLW ), Birk s et al. (120061 hereafter BFG), 
and Lopez-Sepulcre et al. (|2p09). The 12 CO(1-0) line emission 
appears to trace also a N-S collimated flow (see BFG), coincid- 
ing with a thermal radio jet (Heaton & Littl e [19881 GHLW) seen 
also at IR w avelen gths (Dent et al.|1985bl Walth er et a l. 119901 
Fuller et al. I200T1 De Buizer [20061 Zhang et al. I20L31 . It has 
been proposed that the poorly collimated NE-SW outflow and 
the N-S jet could be manifestations of the same bipolar flow un- 
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Fig. 1. Large scale image of the 4.5 fim emission towards the 
star forming region G35.20-0.74N overlaid with a contour map 
of the 350 GHz continuum emission obtained with ALMA. 
The IR image has been obtained applying HiRes deconvolution 
(Velusamy et al. 120081) to the Spitzer/IRAC data. The sub-mm 
map has been corrected for primary beam attenuation. The con- 
tours range from 9 to 203 mJy/beam in steps of 10.8 mJy/beam. 
The dashed circle denotes the primary beam of the ALMA 12-m 
antennas. The ellipse in the bottom right is the ALMA synthe- 
sized beam. 

dergoing precession (Little et al. 1998 1. However, evidence for 
multiple outflows in this region is provided by SiO, H 13 CO + , 
H 13 CN, and H 2 line observations (GHLW; Lee et al. l20T2l . 

A molecular clump elongated perpendicular to the NE-SW 
outflow has been mapped in dense gas tracers (NH3, CS), whose 
emission exhibits a velocity gradient from NW to SE (Little et 
al. 119851 Brebner et al. 1987 >. This was first interpreted as a 
large (~1' or 0.6 pc) disk/toroid rotating about the NE-SW out- 
flow axis, but GHLW, on the basis of their H 13 CO + and H 13 CN 
observations, propose that this is actually a fragmented rotat- 
ing envelope containing multiple young stellar objects (YSOs). 
Indeed, GHLW identify a core at the center of the outflow and 
another core, named G35MM2, offset to the SE. 

2. Observations and results 

G35. 20-0.74 N was observed with ALMA Cycle at 350 GHz 
in May and June 2012, with baselines in the range 36^-00 m, 
providing sensitivity to structures <2 ". The digital correlator 
was configured in 4 spectral windows (with dual polarization) of 
1875 MHz and 3840 channels each (covering the ranges 334.85- 
338.85 GHz and 346.85-350.85 GHz), providing a resolution of 
~0.4 km s -1 . Flux, gain, and bandpass calibrations were obtained 
through observations of Neptune and J 175 1+096. The data were 
calibrated and imaged using CASA. A continuum map was ob- 
tained from line-free channels and subtracted from the data. The 
synthesized beam is 0751 x 0746, P. A. =48°. The rms noise is 
~6 mJy beam -1 for individual line channels, while in the contin- 
uum image it is ~1.8 mJy beam -1 , implying a S/N of only ~100. 
The latter indicates a reduced dynamic range. 

In Fig. [T| we present the map of the 350 GHz continuum 
emission overlaid on an enhanced resolution Spitzer/IRAC im- 
age at 4.5 /zm extracted from the GLIMPSE survey (Benjamin 
et al. 120031) . The sub-millimeter continuum emission is clearly 
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Fig. 2. Same as Fig. [T] with overlaid the 3.6 cm continuum map 
by GHLW (magenta contours). Unlike Fig. Q] the 4.5 fim im- 
age is displayed using linear scaling to emphasize the brightest 
structures outlining the N-S jet. Cores A, B, and G35MM2 are 
indicated by the arrows. 



tracing an elongated structure across the waist of the butterfly 
shaped nebula. In all likelihood, we are detecting the densest part 
of the flattened molecular structure observed on a larger scale by 
Little et al. ( fl"985i l, Brebner et al. ([T987j l, and GHLW. 

Along the elongated structure a chain of at least 5 cores is 
seen (Fig. [2]), lending support to GHLW's idea that one is deal- 
ing with a fragmented structure instead of the smooth disk/toroid 
hypothesized by Little et al. (1985). We stress that the an- 
gular resolution of our maps (~7 times better than previous 
(sub)millimeter observations) reveals that the YSOs powering 
the outflow(s) lie inside cores A and/or B (see Fig. [2]), because 
these two are the only cores located close to the geometrical cen- 
ter of the bipolar nebula. In particular, core B lies along the N-S 
jet traced by the IR and radio emission, and coincides with one 
of the free-free sources detected by GHLW. The hypothesis by 
BFG that core G35MM2 could be driving the NE-SW outflow is 
not convincing, as this core is off-center, lying right at the border 
of the waist outlined by the IR emission. 

Methyl cyanide emission (as well as other typical HMC 
tracers) is clearly detected only towards cores A and B, and 
marginally towards G35MM2. Emission by vibrationally excited 
lines of CH3CN also indicates that cores A and B could be host- 
ing massive stars. Core B coincides also with a compact free- 
free continuum source detected by GHLW at 6 and 3.6 cm, and 
by Codella et al. (120101 ) at 1 .3 cm. This emission could be part 
of the N-S thermal radio jet or might be coming from an H11 
region ionized by an embedded early-type star. We will discuss 
this possibility in Sect. [4] Faint radio emission at a 5.6<r level 
is detected also towards core A, consistent with the presence of 
embedded star(s). 

We now investigate the gas velocity field in the two cores 
by computing the first moment of a CH3CN line, a dense gas 
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Fig. 3. Overlay of the average emission in the CH3CNXI9-I8) 
K-2 line (white contours) and 3.6 cm free-free continuum emis- 
sion (red contours; from GHLW), on a map of the first moment 
(i.e. velocity) of the same CH3CN line. The velocity intervals 
used to calculate the mean emission and first moment are 25.9- 
38.0 km s _1 for core A and 25.9-33.0 km s for core B. The 
dashed line is intended to stress that different colour scales are 
used for the two cores to emphasize the velocity gradients. The 
white contour levels range from 80 to 656 mJy/beam in steps of 
96 mJy/beam. The crosses mark the positions of the 350 GHz 
continuum peaks. The ellipse in the bottom right indicates the 
synthesized beam. 



tracer. Figure |3]plots the result for the CH 3 CN(19-18) K=2 line, 
with overlaid the line emission averaged over the same velocity 
interval used to calculate the first moment. Note that the mean 
velocity of core B (~30 km s _1 ) differs by ~2 km s _1 from that 
of core A (~32 km s -1 ), which also differs by a similar amount 
from that (~34 km s _1 ) obtained with lower angular resolution 
observations of other tracers (see e.g. GHLW). Discrepancies 
of this type are often found in high-mass star-forming regions 
(see, e.g., th e case of AFGL5142 - Estalella et al. [19931 Zhang 
et al. 120071) and are likely related to the fragmentation process 
in molecular clumps. The interesting result is that both cores 
present velocity gradients, which may be due to expansion or 
rotation of the gas. Note that the sense of the gradient in core B 
is the opposite of that measured on larger scales (see Fig 5 of 
Little et al. [19851 and Figs. 8 and 9 of GHLW). However, the ve- 
locity range traced by the extended emission (30 — 41 km s -1 ) dif- 
fers significantly from that sampled by the CH3CN lines (~26— 
38 km s _1 ), suggesting that we are observing cores whose veloc- 
ity field is not tightly related to that of the largest, fragmenting 
structure where they are embedded. In the following, we will 
address this question for core B, which is better resolved and 
associated with the N-S jet/outflow. 



3. The nature of core B: A Keplerian rotating disk? 

An elongated structure with a velocity gradient along its major 
axis, such as that observed in core B, can be interpreted as either 
a collimated bipolar outflow, or a rotating ring/disk seen close 
to edge on. We thus first investigate the possibility that we are 
observing the root of a large-scale bipolar outflow. 



At first glance, both the orientation (see Fig. [2]i and sign 
of the velocity gradient seen in Fig. [3] appear consistent with 
those of the N-S jet imaged by BFG (see the red-shifted emis- 
sion at 40 km s in their Fig. 3). However, the CO emission in 
the northern jet lobe is detected by BFG also at 28 km s _I (see 
their Fig. 3), namely at blue-shifted velocities. Moreover, the di- 
rection of the CH3CN velocity gradient has a negative position 
angle, whereas the jet is slightly inclined to the east. Such a dif- 
ference could be explained if the jet is precessing (as suggested 
by BFG) and lies close to the plane of the sky, because the jet 
direction would change from the small to the large scale. We 
believe that this explanation is not satisfactory, because the ori- 
entation of the CH3CN velocity gradient differs from that of the 
radio jet on the same scale. This can be seen in Fig. [3] where the 
two free-free emission sources are aligned N-S, whereas the ve- 
locity gradient (and the core major axis) are inclined to the west 
(PA^-20°). 

Based on these findings, we consider the alternative possi- 
bility that the velocity gradient is tracing a disk. To analyse the 
kinematics of the gas in better detail, we fitted the CH3CN emis- 
sion in each 0.4 km s velocity channel with a 2-D Gaussian. 
This allows us to derive the peaks of the line emission at differ- 
ent velocities. Then all peaks can be plotted together to obtain a 
picture of the gas velocity field. This is done in Fig. [4^, where 
for each peak we also plot the corresponding 50% contour level 
to give an idea of the size of the gas emitting at that velocity. 

Two considerations are in order. The first is that the peaks 
outline a sort of elliptical pattern, roughly centered on the posi- 
tion of the continuum peak and oriented SSE-NNW. The second 
is that the most blue- and red-shifted peaks tend to converge to- 
wards the position of the sub-mm continuum peak. On this basis, 
one is tempted to hypothesize that the CH3CN emission is trac- 
ing Keplerian-like rotation, as this would explain why both the 
emission at systemic velocities and that at high (blue- and red- 
shifted) velocities are observed towards the rotation center. 

With this in mind, we have fitted the observed velocity pat- 
tern assuming Keplerian rotation about a point-like source. We 
stress that ours is a purely kinematical fit, where only the rota- 
tion velocity in a 2-D disk is calculated, and no estimate of the 
observed line intensity is computed. Also, no attempt is made to 
demonstrate the uniqueness of the model proposed and we envis- 
age the possibility that other interpretations might be possible. 
Our simple model is sufficient to constrain a number of impor- 
tant parameters. The inputs of the model are: the LSR velocity 
of the system (Vs ys ), the central point mass (M.) and its position 
(xo,yo), the inclination angle (ff) of the rotation axis with respect 
to the plane of the sky, and the position angle (tfr) of the projected 
major axis of the disk. 

The best fit was obtained by varying all input parame- 
ters inside reasonable ranges and minimizing the expression 

£ .^observed _ ymodel^ 2 where • indicates the peak m a generic 

velocity channel. For our calculations, we have used the K-2, 
3, and 4 components of the CH3CN(19-18) transition as well 
as other prominent lines of other species, such as CH30H(7i 6- 
61,5) vt=l, CH 3 OH(14 M 3-14o, M ). and HC 3 N(37-36). All of 
these lines reveal a kinematical pattern very similar to that in 
Fig- Ek' demonstrating that such a pattern does not depend on 
the tracer and mirrors a real physical structure. In Fig. 5J> we 
show an overlay of the velocity peaks of all the lines and a 
velocity map of the best-fit model. Clearly, the agreement be- 
tween the computed and observed LSR velocities is remarkable. 
The best-fit parameters are the following: V sys =30.0+0.3 km s _1 , 
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Fig. 4. (a) Peaks of the 
CH 3 CN(19-18) K=2 lines 
emission (solid circles) ob- 
tained with a 2-D Gaussian fit 
channel by channel. For each 
peak also the corresponding 
50% contour level is drawn 
(using the same colour as the 
peak). The colour corresponds 
to the line-of-sight velocity, ac- 
cording to the scale displayed 
to the right. The ellipse in the 
bottom right denotes the syn- 
thesized beam, (b) Comparison 
between the velocity of the 
best-fit Keplerian disk (colour 
map) and the emission peaks 
at different velocities (solid 
circles) obtained for the fol- 
lowing lines: CH 3 CN(19-18) 
£=2,3,4, CH 3 OH(7i 6 -6i 5 ) 
v t =l, CH 3 OH(14 U3 -i4o,i 4 ), 
HC 3 N(37-36). The crosses 
mark the position of the 
continuum peak. 



x =18 h 58 m 13 s .027+0 i .002, yo=01 40'35'.'94±0'.'07, iA=157±4°, 
0=19+1°, M t =18+3M . 

A lower limit for the disk radius is — 2500 AU, obtained 
from the maximum deprojected distance of the peaks from the 
center. One may wonder if such a big structure can undergo 
Keplerian rotation. The mass of the disk can be estimated from 
the continuum emission. The integrated flux density from core 
B at 350 GHz is 0.32 Jy. Assuming a dust absorption coefficient 
0.5 cm 2 g^ (v/230.6 GHz) (Kramer et al. 120031 and a gas-to- 
dust mass ratio of 100, we obtain ~3 M for a dust temperature 
of 100 K. Despite the large uncertainties on the dust opacity and 
temperature, we believe that the mass of core B is significantly 
less than the central mass (~18 M ), which satisfies the condi- 
tion for Keplerian rotation. 

It is worth stressing that our findings are in good agreement 
with the recent study by Ilee et al. (120121 1. Through measure- 
ments of scattered light from G35.20-0.74N, these authors find 
that the CO first overtone bandhead emission at 2.3 jum can be 
fitted with a Keplerian disk rotating about a 17.7 M Q star. 

The distribution of the molecular peaks in Fig. |4j) clearly 
shows that our observations detect only the NE side of the disk. 
We speculate that this could be an opacity effect. If the disk is 
optically thick in the relevant lines, flared, and inclined by 19°, 
only part of the surface heated by the star is visible to the ob- 
server. This creates an asymmetry along the direction of the pro- 
jected disk axis, with line emission being more prominent on the 
side (in our case the NE side) where the disk surface is visible. 
Clearly, radiative transfer calculations are needed to confirm this 
scenario, but we note that in our source the NE part of the disk 
axis should be pointing towards the observer, consistent with the 
orientation of the CO outflow (blue shifted to the NE and red 
shifted to the SW - see GHLW) and the obscuration seen to the 
SW in the IR images (see Fig. 2). 



4. The stellar content of core B: A binary system? 

An issue that is worth discussing is whether an 18 M YSO is 
compatible with the bolometric luminosity of the region. (3 x 
10 4 L ; see Sect.Q]). Depending on the adopted mass-luminosity 



relation, the luminosity expected for an 18 M star ranges from 
2.5 x 10 4 L Q (Diaz-Miller et al. ll998l . to 6.6 x 10 4 L (Martins 
et al. 120051 . This means that the 18 M star should be the main 
contributor to the luminosity of the whole star forming region. 
Such a possibility seems quite unlikely due to the presence of 
multiple cores (see Fig. |5J, one of which is an HMC possibly 
hosting at least another high-mass star (core A). 

A possibility is that one is underestimating the true lumi- 
nosity due to the "flashlight effect", where part of the stellar 
photons are lost through the outflow cavities. According to the 
recent model by Zhang et al. (120131 1. when this effect is taken 
into account, the luminosity obtained assuming isotropic emis- 
sion (3.3 x 10 4 L s ) becomes as large as 7 x 10 4 -2.2 x 10 5 L 
consistent with a single star of ~20-34 M e . 

While the previous explanation is possible, the isotropic es- 
timate appears more robust than a model-dependent value, and 
we thus consider another hypothesis, namely that one is dealing 
with a binary system. In this case, the luminosity is significantly 
reduced with respect to that of a single 18 M object and may 
be as low as ~ 7 x 10 3 L Q for equal members. The latter is much 
less than the bolometric luminosity, thus allowing for a signifi- 
cant contribution by other YSOs. 

The existence of a binary system could also explain why 
the N-S jet associated with core B is not aligned with the disk 
rotation axis. The presence of a companion would in fact be 
sufficient to induce precession of the jet/outflow, as hypothe- 
sized by Shepherd et al. (2000) to explain the observed pre- 
cession of the bipolar outflow from the high-mass protostar 
IRAS 20126+4104. In this scenario, the outflow from core B 
would precess about an axis oriented NE-SW, i.e. along the bi- 
sector of the butterfly shaped IR nebula seen in Fig. Q] The IR 
emission would arise from the cavity opened by the outflow it- 
self during its precession, while the thermal radio jet would trace 
the current direction of the precessing axis. 

The last question we address is the origin of the free-free 
emission from core B (see e.g. Fig. |3j. Could this be tracing 
a (hypercompact) Hn region? According to GHLW, this source 
(n. 7 in their Table 1) has a spectral index >1.3 between 6 
and 3.6 cm, compatible with free-free emission from an opti- 
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cally thick Hn region. Extrapolation of the 3.6 cm flux den- Zhang, Q., Hunter, T. R., Beuther, H., et al. 2007, ApJ, 658, 1152 
sity (0.5 mJy) to 1.3 cm gives >1.9 mJy, in agreement with the Zinnecker, H. & Yorke, H. w. 2007, ARA&A, 45, 481 
3 mJy flux measured, with lower angular resolution, by Codella 
et al. (1201 Oi l. It is hence possible that the emission is partially 
thick at 1.3 cm and the Lyman continuum estimate of 10 45 s _1 
obtained under the optically thin assumption is a lower limit. A 
binary system with a total mass of 18 M & has a Lyman contin- 
uum flux (see Diaz-Miller et al. 1998 ) ranging from 5 x 10 44 s _1 
(for equal masses) to 1.8xl0 47 s _1 . We conclude that we could be 
observing an Hn region ionized by a binary system at the center 
of a Keplerian disk. Should this be confirmed, G35.20-0.74N 
would represent a unique example of radio jet coexisting with 
an Hn region powered by the same YSO(s). Investigating such a 
short-lived transition phase in the evolution of an OB-type star 
could provide us with important clues on the formation process 
of these objects. 
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